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ABSTRACT 

As  Vertical  Line  Arrays  (VLA)  are  used  at  increasingly  higher  frequencies  (3  kHz),  the 
importance  of  the  straightness  of  the  array  increases.  Finite  Element  Analysis  (FEA)  has 
been  used  to  estimate  the  shape  of  typical  VLAs  subject  to  currents  in  shallow  water. 
Two  typical  VLA  configurations  have  been  modelled  and  the  results  presented. 
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Estimating  the  Straightness  of  Vertical  Line 
Arrays  using  Finite  Element  Analysis 


Executive  Summary 

Vertical  Line  Acoustic  Arrays  (VLA)  are  being  used  at  increasingly  higher  frequencies 
(3  kHz).  At  these  frequencies,  the  alignment  of  the  hydrophones  in  die  array  relative  to 
each  other  (straightness)  is  of  great  importance.  One  method  of  straightening  the  array 
is  to  increase  the  tension  applied  to  it.  Finite  Element  Analysis  (FEA)  has  been  used  to 
estimate  the  straightness  of  two  typical  VLAs  subject  to  current  flow  in  shallow  waters. 

The  analysis  shows  that  the  position  of  the  array  can  be  described  as  a  horizontal 
displacement  of  a  reference  point,  a  rotation  or  tilt  about  that  point  and  the  deviation  of 
the  hydrophones  from  the  mean  line  through  the  hydrophone  positions.  In  terms  of 
beamforming,  displacement  should  cause  no  error  except  in  the  very  near  field; 
rotation  causes  some  error  and  deviation  the  greatest  error. 


FEA  demonstrates  that  adding  tension  to  the  VLA  decreases  displacement,  rotation 
and  deviation.  However,  there  is  a  practical  limit  to  the  reduction  obtainable  as  each 
variable  asymptotes  rapidly  towards  a  constant  value  as  the  tension  is  increased. 
Adding  tension  to  the  VLA  is  most  efficient  at  reducing  displacement  and  least 
efficient  at  reducing  deviation. 

The  conclusion  of  the  analysis  is  that  applying  tension  to  a  VLA  is  not  an  effective  or 
practical  method  of  reducing  the  rotation  and  deviation  of  the  array  when 
beamforming  at  frequencies  in  the  3  kHz  band. 
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1.  Introduction 


As  Vertical  Line  Arrays  (VLA)  are  used  at  increasingly  higher  frequencies  (3  kHz),  the 
importance  of  tire  straightness  of  the  array  increases.  Finite  Element  Analysis  (FEA)  has 
been  used  to  estimate  the  shapes  of  typical  VLAs  subject  to  currents  in  shallow  water. 
Two  typical  arrays,  one  bottom  moored  the  other  surface  moored,  have  been  modelled 
and  the  results  presented. 


2.  The  Models 

The  computer  models,  illustrated  in  Figure  2.1,  were  generated  and  analysed  using  the 
ANSYS®  finite  element  analysis  software  [1]  run  on  a  Unix  workstation.  Sample  model 
and  solution  input  files  for  tire  Bottom  Moored  Array  are  presented  in  Appendix  A  and 
the  Surface  Moored  Array  in  Appendix  B.  The  models  are  highly  non-linear  and  both 
stress  stiffening  and  large  deflection  effects  have  to  be  included.  All  external 
components  were  modelled  using  the  PIPE59,  'immersed  pipe  or  cable'  element.  This  is 
the  only  element  that  can  include  the  effects  of  external  loads,  inertia,  buoyancy  and 
hydrodynamic  forces.  The  tube  of  an  array  is  usually  filled  with  an  'incompressible' 
fluid  and  pressurised  to  approximately  1  to  2  atmospheres.  As  this  element  cannot 
model  a  contained,  incompressible  fluid,  an  increased  internal  pressure  of  4 
atmospheres  was  applied  to  the  tube  to  avoid  'tube  collapse'  error  messages  that 
frequently  occur  during  tire  solution  sequence.  All  external  elements  were  assigned  a 
Reynolds  Number  independent,  normal  drag  coefficient  of  1.2.  Internal  braided  cord 
strength  members  were  modelled  using  the  non-linear  UNK10,  'tension  only  spar' 
element.  The  models  were  considered  to  be  moored  in  40  m  depth  of  sea  water  with  a 
density  of  1025  kg/m3. 


Figure  2.1  Schematic  of  Models 
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2.1  Bottom  Moored  Array 

The  first  VLA  modelled  was  bottom  moored  and  comprised  a  float,  float  cable,  array 
and  mooring  cable  fixed  at  the  mudline.  The  complete  VLA  has  a  mass  of  162.3  kg  in 
air  and  has  a  negative  buoyancy  force  of  15.9  N  in  seawater.  The  mechanical  and 
physical  properties  of  the  VLA  are  shown  in  Table  2-1  and  Table  2-2. 


Table  2-1  Mechanical  Properties  -  Bottom  Moored  Array 


Modulus 
of  Elasticity 
[Pa] 

Density 

Poissons 

Ratio 

Internal 

Components 

fkg/m] 

Float 

Float  Cable 

Strength  Member 

Array  Tube 

Mooring  Cable 

68.26x10’ 

20.0x10’ 

20.-0x10’ 

18.4  x  106 
53.3x10’ 

1025 

763 

1067 

1320 

5385 

0.33 

0.29 

0.33 

0.33 

0.29 

0.8186 

Table  2-2  Physical  properties  -  Bottom  Moored  Array 


Length 

Diameter 

Wall 

Thickness 

Cross 

Sectional 

Area 

Element 

Length 

Number 

of 

Elements 

Crossflow 

Drag 

Coefficient 

m 

[mm] 

[mm] 

fm2] 

M 

Float 

Float  Cable 

Strength  member  #1 
Strength  member  #2 
Array  Tube 

Mooring  Cable 

1 

4 

24 

24 

24 

8 

400 

1Z7 

40 

8 

3 

4.28  x  Id6 
4.28  x  10-* 

1 

1 

1 

1 

4 

24 

24 

24 

8 

1.2 

1.2 

1.2 

1.2 

2.1.1  Float 

To  enable  the  tension  in  the  array  to  be  easily  changed,  the  float  was  modelled  as  a 
rigid  cylinder  of  neutral  buoyancy,  and  a  vertical,  upward  force  was  applied  at  the  top 
of  the  float  to  simulate  its  buoyancy  force. 


2.1.2  Float  Cable 

The  float  cable  was  modelled  as  a  polypropylene  rope. 

2.1.3  Acoustic  Array 

The  acoustic  array  comprised  an  outer  cylindrical,  hollow  tube,  modelled  using  the 
cable  formulation  (no  bending  stiffness),  containing  two  chord-like  strength  members. 
The  strength  members  were  constrained  in  position  in  the  tube  at  1  m  intervals  as 
illustrated  in  Figure  2.2.  The  mass  of  the  internal  components  and  fluid  fill  was  set  to 
make  the  array  neutrally  buoyant. 
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Figure  2.2  Schematic  of  Array  Model 
2.1.4  Mooring  Cable 

The  mooring  cable  was  modelled  as  a  stranded  steel  cable.  It  was  anchored  at  the 
bottom. 

2.2  Surface  Moored  Array 

The  second  VLA  modelled  was  surface  moored  and  comprised  a  mooring  rope, 
electronics  can,  array  and  a  weight.  The  complete  VLA  has  a  mass  of  259  kg  in  air  and 
has  a  positive  buoyancy  force  of  2.9  N  in  sea  water.  The  mechanical  and  physical 
properties  of  the  VLA  are  shown  in  Table  2-3  and  Table  2-4. 

2.2.1  Mooring  Rope 

The  mooring  rope  was  modelled  as  a  polypropylene  rope.  It  was  anchored  at  the 
surface. 

2.2.2  Electronics  Can 

The  electronics  can  was  modelled  as  a  rigid  cylinder  of  neutral  buoyancy. 

2.2.3  Acoustic  Array 

The  acoustic  array  comprised  an  outer  cylindrical,  hollow  tube  containing  two  strength 
members,  fluid  fill  and  internal  components.  The  strength  members  were  constrained 
in  position  in  the  tube  at  approximately  1  m  intervals.  The  mass  of  the  internal 
components  and  fluid  fill  was  set  to  make  the  array  neutrally  buoyant. 
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Table  2-3  Mechanical  Properties  Surface  Moored  Array 


Modulus 

[Pa] 

Density 

[kg/m33 

Poissons 

Ratio 

Internal 

Components 

[kg/m] 

Mooring  Rope 
Electronics  Can 
Strength  Member 

Array  Tube 

Weight 

20.0  x  10* 

68.3  x  109 

20.0  x  10» 

18.4  x  106 

68.3  x  109 

763 

1025 

1067 

1320 

1025 

0.29 

0.33 

0.33 

0.33 

0.33 

0.8186 

Table  2-4  Physical  Properties  -  Surface  Moored  Array 


Length 

m 

Diameter 

[mm] 

Wall 

Thickness 

[mm] 

Cross 

Sectional 

Area 

M 

Element 

Length 

[m] 

Number 

of 

Elements 

Crossflow 

Drag 

Coefficient 

Mooring  Rope 

9 

12.7 

1 

n  ^ 

9 

1 

1.2 

1.2 

Electronics  Can 

0.5 

220 

U«J 

Strength  member  #1 

28.5 

4.28  x  10* 

A  OQ  *1 

31 

31 

Strength  member  #2 

28.5 

4.2o  X  It/* 

31 

1.2 

Array  Tube 

28.5 

40 

3 

1 

1 

1.2 

Weight 

1.0 

500 

2.2.4  Weight 


To  enable  die  tension  in  the  array  to  be  easily  changed,  the  weight  was  modelled  as  a 
rigid  cylinder  of  neutral  buoyancy,  and  a  vertical,  downward  force  was  applied  at  the 
bottom  of  the  weight  to  simulate  its  inertia  (gravitational)  force. 


3.  Loads 


3.1  Inertia  Forces 

Both  models  were  loaded  by  inertia  forces  generated  by  the  software  from  the  masses 
and  acceleration  due  to  gravity  (9.81  m/s2)  acting  vertically  down. 

3.2  Buoyancy  Forces 

The  software  generated  the  appropriate  buoyancy  forces  to  be  applied  to  each  model 
using  the  geometric  and  density  data. 

3.3  Hydrodynamic  Forces 

Both  models  had  hydrodynamic  forces  generated  by  the  software  using  die  geometry, 
drag  coefficients  and  the  current  velocity.  Two  current  profiles  were  used. 
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3.3.1  Constant  Current 

The  applied  current  was  considered  constant  throughout  the  water  depth  at  0.6  knots 
(0.31  m/s). 

3.3.2  Profiled  Current 

The  applied  current  was  shaped  with  water  depth  as  shown  in  Figure  3.1. 


Figure  3.1  Shaped  Current  profile 

3.4  Applied  Tension 

An  applied  force  tensioned  the  VLA.  This  force  was  varied  to  illustrate  the  variation  of 
the  straightness  of  the  array  with  the  applied  tension.  The  force  ranged  from  100  N  to 
1000  N  in  steps  of  100  N  in  most  runs. 

3.4.1  Bottom  Moored  Array 

The  force  was  applied  vertically  upward  at  the  top  of  the  float.  This  simulated  a  float 
with  fixed  geometry,  hence  constant  drag  forces,  and  variable  buoyancy. 

3.4.2  Surface  Moored  Array 

The  force  was  applied  vertically  downward  at  the  bottom  of  the  weight.  This  simulated 
a  weight  with  fixed  geometry,  hence  constant  drag  forces,  and  variable  mass. 
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4.  Results 


The  FEA  models  were  run  using  ANSYS  finite  element  analysis  program.  Results  were 
either  plotted  directly  in  ANSYS  or  transferred  to  EXCEL  for  further  analysis  and 
plotting.  Table  4-1  summarises  the  parameters  used  in  each  analysis. 


Table  4-1  Summary  of  Analyses 


Analysis 

Folder 

Model 

File 

Solution 

File(s) 

Tension  Range 
[N] 

Min.  Max. 

Sea 

State 

Min 

Current  R 
[knofc 
Max 

a 

spreadl 

solOl+02 

1000 

0 

0.6 

b 

spreadl 

sol03 

20 

1000 

0 

0.6 

c 

spreadl 

sol03 

100 

1000 

0 

0.6 

d 

spreadl 

sol04 

100 

1000 

0 

0.6 

e 

spreadl 

solOS 

100 

1000 

0 

0.6 

f 

spreadla 

so!04 

100 

1000 

0 

0.6 

g 

spreadl 

sol06 

100 

1000 

3 

‘  0.6 

h 

spreadl 

sol09 

100 

1000 

3 

0.6 

i 

spreadl 

sol07 

500 

0 

0.1 

1.0 

a 

spread2 

solution.  01 

100 

1000 

0 

0.6 

Profile 


Comment 


Constant 

Constant 

Constant 

Constant 

Shaped 

Constant 

Constant 

Constant 

Constant 

Constant 


Deleted 

Deleted 

Deleted 

Bottom  Moored 

1  Strength  Member 
Deleted 


Surface  Moored 


4.1  Summary  of  Results 

The  deflection  of  the  acoustic  array  can  be  considered  as  having  three  components: 

a)  a  lateral  displacement  of  a  reference  point,  the  node  representing  the 
hydrophone  nearest  the  mooring,  on  the  acoustic  array 

b)  a  rotation  about  the  reference  point  of  the  mean  line  through  the  acoustic 
array  from  the  vertical  and 

c)  a  horizontal  deviation  of  the  nodes  representing  the  hydrophones  from  the 
mean  line. 

Folders  d,  e  and  f  were  analysed  to  determine  the  mean  line  through  the  nodes  in  the 
acoustic  array  at  1000,  500  and  100  N.  The  mean  line  was  obtained  by  a  linear  'least 
squares  fit'  through  the  deflected  position  of  the  nodes.  The  rotation  of  this  mean  line 
to  the  vertical  and  the  horizontal  deviation  of  the  nodes  from  this  line  were  calculated 
The  total  deviation  is  the  sum  of  the  maximum  positive  and  negative  deviations 
obtained.  The  results  are  summarised  in  Table  4-2  and  Table  4-3. 


Examination  of  these  tables  shows  that  the  maximum  rotation  and  deviation  occurs 
with  a  constant  current  profile.  The  shaped  current  profile  causes  lower  hydrodynamic 
loads  due  to  the  decreased  mass  flow  past  die  VLA,  although  the  maximum  current 
speed  was  the  same  in  both  instances. 
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Table  4-2  Results  Summary  -  Bottom  Moored  Array 


Run 


Sea  State 


Maximum  Current 


Speed 

[knots] 


Type 


'ension 

Rotation 

Maximum  Deviation 

-ve 

+ve 

Total 

IN] 

[degrees] 

[m] 

[m] 

[m] 

1000 

3.2 

-0.121 

0.061 

0.182 

500 

6.1 

-0.231 

0.117 

0.348 

100 

19.4 

-0.772 

0.393 

1.165 

1000 

1.7 

-0.106 

0.053 

0.159 

500 

3.2 

-0.206 

0.102 

0.308 

100 

11.8 

-0.766 

0.372 

1.138 

1000 

3.2 

-0.121 

0.061 

0.182 

500 

6.0 

-0.231 

0.117 

0.348 

100 

19.4 

-0.768 

0.391 

1.159 

spreadld 

spreadle 

spreadlf 


0.6 


0.6 


0.6 


Constant 


Shaped 


Constant 


Table  4-3  Results  Summary  -  Surface  Moored  Array 


Run 

Sea  State 

Maximum  Current 

Tension 

Rotation 

Maximum  Deviation 

Speed 

[knots] 

Type 

[N] 

[degrees] 

-ve 

M 

+ve 

[m] 

Total 

[m] 

spread2a 

0 

0.6 

Constant 

1000 

500 

3.4 

6.5 

-0.105 

-0.199 

0.054 

0.102 

0.159 

0.301 

100 

19.8 

-0.622 

0.322 

0.944 

These  tables  also  show  that  although  the  angle  of  rotation  is  greater  for  the  surface 
moored  array  than  for  the  bottom  moored  one  the  deviations  are  less  for  the  former 
mooring  position.  The  differences  in  the  deviations  are  quite  small  and  are  probably 
caused  by  the  geometry  changes  between  the  two  systems. 


4.2  Comparison  of  Bottom  and  Surfaced  Moored  Arrays 

To  compare  the  effect  of  the  mooring  position  runs  spreadld  and  spread2a  are 
considered.  The  results  for  each  type  of  mooring  are  presented  as  five  graphs: 

a)  Horizontal  Displacement, 

b)  Vertical  Displacement, 

c)  Mooring  Forces 

d)  Acoustic  Array  Position, 

e)  Horizontal  Deviation. 

The  latter  two  graphs  (d  and  e)  should  contain  a  family  of  curves,  one  curve  for  each 
tension  considered.  Only  1000,  500  and  100  N  tensions  have  been  computed  and,  as 
these  curves  all  show  similar  trends,  only  the  curves  for  a  tension  of  500  N  are 
presented  here. 

In  discussing  the  graphs  of  horizontal  displacement  alternative  definitions  of  rotation 
and  deviation  are  used.  They  are  less  precise  than  the  original  definition  but  are 
simpler  and  allow  an  appreciation  of  die  deviation  and  rotation  to  be  gleaned  directly 
from  the  horizontal  displacement  graphs. 
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a)  the  rotation  is  determined  from  a  straight  line  through  the  nodes 
representing  the  upper  and  lower  hydrophone  position  of  the  acoustic 
array; 

b)  the  median  is  the  average  horizontal  displacement  of  the  upper  and  lower 
nodes; 

c)  the  deviation  is  tire  horizontal  distance  of  the  centre  node  of  the  acoustic 
array  from  the  median,  this  would  be  zero  for  a  straight  array. 

Using  this  definition,  the  greater  the  separation  of  the  'top'  and  'bottom'  horizontal 
displacement  curves,  the  greater  the  rotation  and  the  greater  the  separation  of  the 
'centre'  curve  from  the  median,  the  greater  the  deviation.  The  rotation  and  deviation 
calculated  by  this  method  are  similar  to,  the  rotation  and  total  deviation  calculated  by 
the  method  described  in  paragraph  4.1  above. 

4.2.1  Bottom  Moored  Array 


Figure  4.1  Horizontal  Displacement  -  Bottom  Moored  Array 

Figure  4.1  shows  the  variation  of  horizontal  displacement  at  the  top,  middle  and 
bottom  of  the  bottom  moored,  acoustic  array  with  tension.  The  median  line  is  included 
for  reference.  As  expected,  the  top  of  the  acoustic  array,  which  is  less  restrained  by  the 
bottom  mooring,  moves  more  than  the  bottom. 
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Comparing  the  top  and  bottom  curves  it  can  be  seen  that  increasing  tension  decreases 
the  separation  and  hence  the  rotation.  It  is  significant  that  the  reduction  in  this 
separation  is  not  as  pronounced  at  tensions  greater  than  500  N 

Comparing  the  median  and  the  centre  curve,  it  can  be  seen  that  increasing  tension 
reduces  the  separation  and  hence  the  deviation.  Again,  it  is  significant  that  this 
deviation  is  sensibly  constant  at  tensions  greater  than  300  N. 

The  form  of  these  curves  indicates  that  there  is  an  economic  limit  to  the  tension  applied 
in  terms  of  both  array  rotation  and  deviation.  In  this  instance  increasing  the  tension 
beyond  300  N  is  uneconomic  and  would  result  in  a  bulky,  expensive  and  hard  to 
handle  system. 


Figure  4.2  Vertical  Displacement  -  Bottom  Moored  Array 


Figure  4.2  shows  the  variation  with  tension  of  vertical  displacement  of  the  top,  centre 
and  bottom  of  the  acoustic  array.  The  vertical  displacement  is  relatively  small 
compared  with  the  horizontal  displacement.  The  major  component  of  this 
displacement  is  caused  by  the  inclination  of  the  VLA  away  from  the  vertical,  thus,  as 
the  tension  is  increased  the  angle  is  reduced  and  die  array  moves  towards  the  surface. 
The  small  positive  displacement  at  the  higher  tensions  is  caused  by  stretch  of  the  array 
components. 
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Figure  4.3  Mooring  Forces  -  Bottom  Moored  Array 

Figure  4.3  shows  the  variation  of  the  horizontal  and  vertical  components  of  the 
mooring  load  with  tension.  As  expected,  the  vertical  component  closely  follows  the 
tension  applied  to  the  float.  The  horizontal  component,  due  mainly  to  the  crossflow 
drag,  is  substantially  constant  except  at  very  small  tensions  when  it  reduces  due  the 
effect  of  increasing  array  rotation. 

Figure  4.4  shows  the  original  and  deflected  position  of  the  acoustic  array  nodes  due  to 
the  current  flow  when  subjected  to  a  tension  of  500  N.  Also  shown  are  the  best  fit 
straight  line  and  its  rotation  from  the  vertical. 

Figure  4.5  shows  the  horizontal  position  error  of  the  nodes  from  the  best  fit,  straight 
line.  The  maximum  negative  error  occurs  at  the  top  of  the  acoustic  array.  The 
maximum  positive  error  occurs  near  die  centre  of  the  acoustic  array. 

4.2.2  Surface  Moored  Array 

Figure  4.6  shows  the  horizontal  displacement  at  the  top,  centre  and  bottom  of  the 
surface  moored,  acoustic  array.  The  curves  are  similar  in  shape  to  the  bottom  moored 
array  except,  as  expected,  the  bottom  of  the  acoustic  array  moves  more  than  the  top 
which  is  restrained  by  the  surface  mooring. 
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Figure  4.4  Acoustic  Array  Position,  500  N  Tension  -  Bottom  Moored  array 
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Figure  4.5  Acoustic  Array  Horizontal  Deviation,  500  N  Tension  -  Bottom  Moored  Array 
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4.6  Horizontal  Displacement  -  Surface  Moored  Array 


Figure  4. 7  Vertical  Displacement  -  Surface  Moored  array 
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Figure  4.7  shows  the  variation  with  tension  of  vertical  displacement  of  the  top,  centre 
and  bottom  of  the  surface  moored,  acoustic  array.  As  with  the  bottom  moored  array, 
(Figure  4.2)  the  vertical  displacement  is  relatively  small  compared  with  the  horizontal 
displacement.  The  major  component  of  this  displacement  is  caused  by  the  inclination  of 
the  VLA  away  from  the  vertical,  thus,  as  the  tension  is  increased  the  angle  is  reduced 
and  the  array  moves  towards  the  bottom  The  small  negative  displacement  at  the 
higher  tensions  is  caused  by  stretch  of  the  array  components. 


Figure  4.8  Mooring  Forces  -  Surface  Moored  Array 

Figure  4.8  shows  the  variation  of  the  horizontal  and  vertical  components  of  the 
mooring  load  with  tension.  As  expected,  the  vertical  component  closely  follows  the 
tension  applied  to  the  float.  The  horizontal  component,  due  mainly  to  the  crossflow 
drag,  is  substantially  constant  except  at  very  small  tensions  when  it  reduces  due  the 
effect  of  increasing  array  angle.  The  offset  data  points  between  500  and  700  N  tension 
appear  to  be  caused  by  a  glitch  in  the  analysis  program  but  this  has  not  been 
confirmed. 

Figure  4.9  shows  the  original  and  deflected  position  of  the  acoustic  array  nodes  due  to 
the  current  flow  when  subjected  to  a  tension  of  500  N.  Also  shown  are  the  'best  fit' 
straight  line  and  its  rotation  from  tire  vertical. 

Figure  4.10  shows  the  horizontal  deviation  of  the  nodes  from  the  straight  line  fit.  The 
maximum  negative  deviation  occurs  at  the  bottom  of  the  acoustic  array.  The  maximum 
positive  deviation  occurs  near  the  centre  of  the  acoustic  array. 
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ure  4.9  Acoustic  Array  Position,  500  N  Tension 
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Horizontal  Deviation  [m] 


Deviation  -  Spt2a-05.xls 


Figure  4.10  Acoustic  Array  Horizontal  Deviation,  500  N  -  Surface  Moored  Array 
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4.3  Reduced  Number  of  Strength  members 

An  analysis  was  carried  out  with  a  bottom  moored  array  with  one  strength  member 
instead  of  two.  Comparing  runs  spreadld  and  spreadlf  in  Table  4-2,  no  appreciable 
change  in  the  results  with  a  single  strength  member,  compared  with  the  two  strength 
member  analysis,  was  observed. 

4.4  Current  Speed 

An  analysis  was  performed,  spreadli,  with  the  bottom  moored  array  to  investigate  the 
effects  of  varying  the  current  speed  from  0  to  1.0  knots.  A  constant  tension  of  500  N 
was  applied. 


Figure  4.11  Horizontal  Displacement,  500  N  Tension  -  Bottom  Moored  Array 


Figure  4.11  shows  the  horizontal  displacement  of  the  top,  centre  and  bottom  of  the 
bottom  moored,  acoustic  array  as  a  function  of  the  current.  The  curves  almost  follow 
the  well-known  square  law  (actually  deflection  =  const  x  current'  9'' ).  It  can  also  be  seen 
that  die  curve  for  die  bodom  and  top  diverge  indicating  an  increasing  angle  to  the 
vertical,  and,  since  the  centre  curve  diverges  from  the  median,  the  curvature  also 
increases  with  current  speed.  Thus,  for  both  types  of  mooring,  it  can  be  assumed  that 
small  increases  in  the  current  speed  will  cause  relatively  large  increases  in  array  angle 
and  curvature. 


15 


DSTO-TN-0197 


4.5  Sea  State 

An  analysis,  spreadlh,  was  performed  to  investigate  the  effect  of  the  sea  state  on  the 
shape  of  the  VLA.  The  model  used  was  the  bottom  moored  array,  sea  state  3  (wave 
amplitude  1.77  m,  wave  period  4.6  s),  a  constant  current  profile  of  0.6  knots  and  the 
applied  tension  was  1000  to  100  N  in  100  N  increments. 

Figure  4 12,  a  reproduction  of  an  output  figure  from  the  ANSYS  software,  shows  the 
horizontal  displacement  with  time.  Table  4-4  provides  a  reference  between  the  time 
axis  and  the  loading  condition.  This  figure  illustrates  the  method  used  to  arrive  at  an 
initial  condition  before  varying  die  tension  loading. 

0  - 100  s  The  gravity,  buoyancy  forces  and  the  tension  load  of  1000  N  is 
gradually  applied  to  the  VLA  over  a  time  of  100  s.  This  provides 
a  steady  state  position  of  the  array  under  'static'  loading.  As  this 
load  is  vertical,  no  horizontal  deflections  are  produced. 

100  -  200  s  The  current  of  0.6  knots  is  applied  as  a  step  input  and  the 
transient  effects  allowed  to  stabilise  over  the  100  s  interval. 

200  -  700  s  A  wave  train  is  applied  to  the  system  and  the  transient  effects 
allowed  to  stabilise  over  the  500  s  interval. 

700  - 1600  s  The  tension  load  is  gradually  reduced  by  100  N  every  100  s. 


Table  4-4  Solution  Pointsjor  Sea  State  Analysis 


Time 

Tension 

[N] 

Current 

[knots] 

0 

0 

0 

100 

1000 

0 

200 

1000 

0.6 

700 

1000 

0.6 

800 

900 

0.6 

900 

800 

0.6 

1000 

700 

0.6 

1100 

600 

0.6 

1200 

500 

0.6 

1300 

400 

0.6 

1400 

300 

0.6 

1500 

200 

0.6 

1600 

100 

0.6 

Sea  State 

0 

0 

0 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


It  can  be  seen  from  Figure  4.12  that  the  introduction  of  the  Sea  State  3  causes  an 
increased  displacement  with  some  increase  in  rotation  and  divergence.  When  the 
reduction  in  tension  is  commenced  at  800  s,  both  the  rotation  and  the  divergence 
increase  with  time  (decreasing  tension).  Thus,  an  increase  in  Sea  State  will  result  in 
greater  angles  and  divergence  than  those  at  sea  state  zero. 
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The  jitter  in  the  curves  is  the  result  of  the  wave  action  and  is  basically  at  the  frequency 
of  the  wave. 


Figure  4.12  Horizontal  Displacement,  Sea  State  3,  Current  0.6  knot  -  Bottom  Moored  Array 
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5.  Conclusions 

Applying  tension  to  limit  the  deflection  of  the  acoustic  section  of  a  VLA  caused  by 
current  flow  is  not  an  effective  method  to  employ  particularly  when  the  required 
deflections  are  small. 

Applying  tension  to  the  array  is  more  effective  in  limiting  the  lateral  displacement  and 
the  rotation,  rather  than  the  deviation.  Referring  to  the  arrays  of  this  report,  we 
conclude  that  500  to  600  N  is  the  practical  limit  for  restricting  the  lateral  displacement 
and  rotation  whereas  300  to  400  N  can  be  considered  the  limit  for  restricting  the 
deviation  (Figure  4.1  and  Figure  4.6). 

These  analyses  have  also  shown; 

a)  there  is  little  difference  in  the  rotation  or  deviation  between  bottom  or 
surface  moored  VLAs; 

b)  for  a  given  tension,  increasing  the  stress  in  the  strength  member(s)  has  no 
effect  on  the  angle  or  deviation; 

c)  assuming  a  fixed  maximum  current  speed,  a  constant  current  profile  will 
always  cause  larger  deflections  than  a  shaped  profile; 

d)  the  relationship  between  deflection  of  the  array  and  current  speed  is 
approximately  a  square  law;  and 

e)  increasing  sea  state  increases  the  displacement,  rotation  and  deviation  of 
the  array.  - 

Methods  that  should  be  considered  for  decreasing  the  angle  and  deviation  of  the 

acoustic  array  are  to:  .  ,  ,  . 

a)  reduce  the  drag  coefficients  of  all  components,  particularly  the  acoustic 

array; 

b)  reduce  the  diameter  of  the  acoustic  array; 

c)  increase  the  bending  stiffness  of  the  acoustic  array;  and 

d)  optimise  the  shape  of  the  float  or  weight  to  reduce  its  drag. 
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Bottom  Moored  Array 
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A.1  Model  Input  File 


!  Model  Designation  -  spreadl 
/TITLE , SPREAD  VLA  Analysis 

i  * 

1  This  model  is  a  24  m  x  40  mm  dia  Vertical  Line  Array 
!  suspended  in  40  m  of  water  in  currents  up  to  0.6  knots. 

!  The  Tension  on  the  VLA  is  controlled  by  an  applied  force  at  the  float 

i  * 


!  Treats  the  tube  of  the  array  as 
I  * 

KEYW,  PR__SET ,  1 
KEYW, PR_STRUC , 1 
/PMETH, OFF 
/PREP7 

i  * 
i  * 

{  Set  Parameters 
rhow  =  1025 
depth  «  40 

i  * 

!  * 

!  Cable  Elements 
ET, 1, PIPE59 
KEYOPT, 1,1,1 
KEYOPT ,1,6,2 
R,l, 0.008,  ,1.2 


RMORE,  ,  , 

o 

o 

H 

UIMP, 1, EX  , 

,  , 5 . 33el0 , 

UIMP , 1 , DENS , 

,  ,5385," 

UIMP , 1 , ALPX , 

,  , 1 . 51e- 05 , 

UIMP , 1 , NUXY , 
TB, WATER, 1 

,  ,0.29 

TBMODIF, 1 , 3 

, depth 

TBMODIF, 1 , 4 

i  * 

ET, 4 , PIPE59 
KEYOPT, 4, 1,1 
KEYOPT, 4, 6, 2 

,  rhow 

R, 4, 0.0127,  , 

1.2 

RMORE,  ,  ,  ,0 

.01 

UIMP, 4, EX  , 

,  ,  20e9 

UIMP, 4, DENS, 

,  ,763 

UIMP, 4, ALPX, 

,  ,1.51e-05 

UIMP, 4, NUXY, 

,  ,0.29 

TBCOPY, WATER, 

i  * 

!  VLA  Elements 
ET, 3 , LINK10 
KEYOPT, 3, 2,2 
KEYOPT, 3, 3, 0 
R, 3,4 .28e-6 

1 , 4 

UIMP, 3, EX  , 

,  , 20e9 

UIMP, 3, DENS, 

,  ,1067, 

UIMP, 3, ALPX, 

,  , 1 . 51e-05 

UIMP,  3, NUXY, 

,  ,0.33 

a  cable. 


!  Water  Density  1025  kg/cu  m 
■  Water  Depth  40  m 


1  Mooring  Cable 
1  Cable  option 

!  Print  member  F  and  M  in  elem  coords 
!  8  mm  Cable,  Cd=1.2 
!  Ct=0.01 


!  Water  Depth 
!  Water  Density 

*  Float  Cable 
!  Cable  option 

!  Print  member  F  and  M  in  elem  coords 
!  12.7  mm  Cable,  Cd=1.2 
!  Ct=0.01 


!  Strength  Member 
!  Small  Stiffnesses  assigned 
I  Tension  Only  (Cable)  Option 
!  Effective  Area  4.28e-6  sq  m 
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ET,2,PIPE59 
KEYOPT, 2 ,1/1 
KEYOPT ,2,6,2 
R,  2, 0.04, 0.003, 1.2 

RMORE,0 .8186,  ,  ,0.01 

UIMP, 2 , EX  ,  ,  , 18 . 4e6 , 
UIMP , 2 , DENS ,  ,  ,1320, 

UIMP, 2 , ALPX,  ,  , 1 . 51e-05 , 
UIMP , 2 , NUXY ,  ,  ,0.33 

TBCOPY, WATER, 1 , 2 

j  * 

I  Float  Elements 
ET, 5 , PIPE59 
R,  5 , 0 . 4 ,  ,1.2 
RMORE,  ,  ,  ,0.1 
UIMP, 5, EX  ,  ,  ,68.26e9 

UIMP , 5 , DENS ,  ,  , rhow 

UIMP, 5, ALPX,  ,  , 1 . 51e-05 , 

UIMP, 5, NUXY,  ,  ,.33 

TBCOPY , WATER ,1,5 

i  * 

!  Define  Nodes 
N, 1,0, 0,-40 
N, 9, 0,0, -32 
Fill, 1,9 
N,  33 , 0, 0, -8 
Fill, 9, 33 
N,  37, 0,0, -4 
Fill, 33, 37 
N, 38, 0,0, -3 

i  * 

!  Define  Elements 
!  Mooring  Cable 
TYPE , 1 , 

MAT, 1 , 

REAL,  1, 

E,  1, 2 

EGEN, 8,1, -1 

i  * 

!  VLA 

TYPE, 2 
MAT,  2 
REAL, 2, 

E, 9 , 10 
TYPE, 3 
MAT,  3 
REAL, 3 
E, 9 , 10 
E, 9 , 10 

EGEN, 24,1, -3 

1  * 

I  Float  Cable 
TYPE, 4 
MAT,  4 
REAL,  4 
E ,  33 , 34 

EGEN, 4,1, -1 

i  * 

!  Float 

TYPE,  5 
MAT,  5 
REAL,  5 
E,  37, 38 
r  * 

!  * 

Pressure=101 . 3E3 
ESEL, S,MAT,  ,2 
CM, ARRAY, ELEM 

SFE , ALL , 1 , PRES ,  , Pressure*4 


!  Tube 

!  Cable  option 

!  Print  member  F  and  M  in  elem  coords 
!  40  mm  x  3  mm  Wall  Tube,  Cd=1.2 
!  Internals  0.8186  kg/m,  Ct=0.01 


!  Float 
!  Ct=0.01 


!  Anchor 

!  Bottom  of  Array 
!  Top  of  Array 
!  Bottom  of  Float 
!  Top  of  Float 


!  Tube 


!  Strength  Member  #1 
!  Strength  Member  #2 


!  Atmospheric  Pressure  [PA] 

!  Select  Array  Elements  Material  2 

!  Form  element  component  ARRAY 

I  Apply  Internal  Pressure  To  ARRAY  Elements 
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ESEL,all  1  Reselect  all  elements 

!  Apply  Constraints 

D, ALL,  UY,0  !  Fix  all  nodes  in  UY 

Dt  1,  UX,0,  ,  ,  ,UZ  “  !  Fix  anchor  node  1  in  UX  and  UZ 

t  * 

!  Load  with  Gravity 
ACEL, 0,0, 9.81 

i  * 
i  * 

/DSCALE, 1 , 1  !  Deflection  Scaling  -  True  Geometry 

/VIEW, 1,  ,-l 

/PNUM, NODE, 1 

/NUMBER, 0 

/PBC,  U,  ,1 

/PBC,  ROT,  ,1 

/PBC, ACEL,  ,1 

EPLO 

SAVE 

i  * 
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A.2  Solution  Input  File 


!  spread. sol 04 
/SOLU 

!  Model  Designation  spreadl 
I  Parameters 
tension=1000 
dtime=100 
current =0 


!  * 

ANTYPE, 4 
TRNOPT, FULL 
LUMPM , 0 

EQSLV, FRONT, le-08 , 0 
AUTOTS , OFF 
NLGEOM, ON 
SSTIF, ON 

nropt,auto 

TIMINT, OFF 
KBC,  0 

i  ★ 

OUTPR,  BASIC,  LAST 
OUTRES ,  ALL ,  ALL 

i  * 

TIME,  dtime 
NSUBST , 1 
NCNV, 2 

F, 38, FZ, tension 
/TITLE , SPREAD  VLA  Analysis 
SOLVE 
1  * 


!  Transient  Effects  OFF 
!  Ramped  Loads 


!  Apply  Tension  to  top  of  cable 
Tension  %tension%  N,  Current  %current%  knots 
I  Solve  Loadstep  #1 


I  * 

current  =  0.6 

I  * 

TB, WATER, 1 

TBMODIF, 2,1, -depth 
TBMODIF, 2,2, current*0 . 5144 
TBCOPY, WATER, 1 , 2 
TBCOPY, WATER, 1 , 4 
TBCOPY, WATER, 1 , 5 

i  * 

AUTOTS, ON 
LNSRCH, ON 
PRED , ON ,  , ON 

TIMINT, ON, STRUC 
DELTIM,0.01, 0,0,0 
NEQIT, 100 
OUTRES , ALL , LAST 
CNVTOL , F ,  ,0.01,2,1 

i  * 

*DO,I, 1,10,1 

dtime=dtime+100 
F, 38, FZ, tension 
TIME, dtime 

/TITLE, SPREAD  VLA  Analysis  - 
SOLVE 

tens ion=tens ion- 1 0 0 


!  0.6  knots 
1  For  Mooring  Cable 


!  For  Tube 
!  For  Float  Cable 
!  For  Float 


!  Turn  ON  Structural  Inertia  Effects 
l  Minimum  Timestep  0.01  s 
1  Maximum  of  100  Equilibrium  Operations 
l  Results  at  end  of  Loadstep 
!  TOLER  =1%,  MINREF  *  1  N 


I  Increment  time 
!  Apply  tension 

Tension  %tension%  N,  Current  % current %  knots 
!  Solve  Loadstep 
!  Decrement  tension 


*ENDDO 

FINISH 
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Surface  Moored  Array 


B.l  Model  Input  File 


!  Model  Designation  -  spread2  model. 01 
/TITLE , Suspended  VIA  Analysis 

i  * 

!  This  model  is  a  Vertical  Line  Array  of  25  Hydrophones  (980  mm  spacing) 
i  Tube  -  28.5  m  Long  x  40  mm  dia 
!  Suspended  in  40  m  of  water 
!  Currents  up  to  0.6  knots. 

!  The  Tension  on  the  VIA  is  controlled  by  an  applied  force  at  the  Weight 


!  The  tube  of  the 

array  is  treated 

j  * 

KEYW, PR_SET, 1 

KEYW, PR  STRUC , 1 

/PMETH, OFF 

/ PREP7 

i  * 

j  * 

!  Set  Parameters 

rhow  =  1025 

i 

depth  -  40 

i 

j  ★ 

ET, 1, PIPE59 

t 

KEYOPT, 1,1,1 

t 

KEYOPT, 1,6,2 

i 

R, 1,0. 0127,  ,1. 

2  1 

RMORE,  ,  ,  ,0.01  ♦ 

UIMP , 1 , EX  ,  , 

,  20e9 

UIMP , 1 , DENS ,  , 

,763 

UIMP, 1, ALPX,  , 

,1.51e-05 

UIMP , 1 , NUXY,  , 

o 

to 

VO 

TB, WATER, 1 

TBMODIF, 1,3, depth  I 

TBMODIF, 1,4, rhow  ! 

!  * 

!  VIA  Elements 

ET , 2 , LINK1 0 

t 

KEYOPT, 2, 2,2 

I 

KEYOPT, 2, 3,0 

{ 

R,  2 , 4 . 28e-6 

i 

UIMP, 2, EX  ,  , 

,  20e9 

UIMP, 2, DENS,  , 

,1067, 

UIMP, 2, ALPX,  , 

,1.51e-05, 

UIMP, 2, NUXY,  , 

,0.33 

ET, 3 , PIPE59 

KEYOPT, 3, 1,1 

KEYOPT ,3,6,2 

R, 3, 0 .04, 0 .003, 

1.2  ! 

RMORE, 0.8186, 

H 

O 

O 

UIMP, 3, EX  ,  , 

, 18 . 4e6 , 

UIMP, 3, DENS,  , 

,1320, 

UIMP, 3, ALPX,  , 

, 1 . 51e-05, 

UIMP, 3, NUXY,  , 

,0.33 

TBCOPY, WATER, 1, 

3 

i  * 

I  Can  Elements 

as  a  cable. 


Water  Density  1025  kg/cu  m 
Water  Depth  40  m 

Mooring  Rope 
Cable  option 

Print  member  F  and  M  in  elem  coords 

12.7  mm  Cable,  Cd=1.2 

Ct=0.01 


Water  Depth 
Water  Density 


Strength  Member 
Small  Stiffnesses  assigned 
Tension  Only  (Cable)  Option 
Effective  Area  4.28e-6  sq  m 


Tube 

Cable  option 

Print  member  F  and  M  in  elem  coords 
40  mm  x  3  mm  Wall  Tube,  Cd=1.2 
Internals  0.8186  kg/m,  Ct=0.01 
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ET , 4 ,  PIPE59 

R, 4 , 0 . 22 ,  ,1.2 
RMORE,  ,  ,  ,0.1 
UIMP,4,EX  ,  ,  , 68  - 26e9 

UIMP , 4 , DENS ,  ,  , rhow 
UIMP , 4 , ALPX ,  ,  ,1.5le-05, 

UIMP, 4 , NUXY,  ,  ,.33 
TBCOPY , WATER ,1,4 

!  * 

!  Weight  Elements 

ET, 5 , PIPE59 
R,  5 , 0 . 5,  ,1.2 

RMORE,  ,  ,  ,0.1 

UIMP, 5, EX  ,  ,  , 68 . 26e9 

UIMP , 5 , DENS ,  ,  , rhow 

UIMP, 5, ALPX,  ,  , 1 . 51e-05, 

UIMP, 5 , NUXY,  ,  ,.33 
TBCOPY, WATER, 1, 5 

i  * 

!  Define  Nodes 
N, 101, 0,0,  0 

N, 110, 0,0,  -9 
Fill, 101,110 
N, 111, 0,0,  -9.5 
N, 115, 0,0, -13 
FILL, 111, 115 
N,  1,0, 0,-14 
N,  25,0,0,-36.5 
Fill, 1,25 
N,201,0, 0,-37 
N,202,0, 0,-38 
N, 203,0, 0,-39 

!  * 

!  Define  Elements 

!  Mooring  Rope 
TYPE , 1 , 

MAT, 1 , 

REAL, 1 , 

E, 101, 102 
EGEN, 9,1, -1 

!  Can 

TYPE, 4 
MAT,  4 
REAL,  4 

E, 110, 111 

i  * 

I  Upper  VLA  -  Can  to  D1 
TYPE , 3 
MAT,  3 
REAL, 3, 

E, 111, 112 
TYPE, 2 
MAT,  2 
REAL,  2 

E, 111, 112 
E, 111, 112 

EGEN, 4, 1, -3 

j  * 

!  Upper  VLA  -  D1  to  HI 
TYPE, 3 
MAT,  3 
REAL, 3, 

E, 115,1 
TYPE, 2 
MAT,  2 
REAL, 2 
E, 115,1 
E, 115,1 


I  Can  500  mm  x  22  mm  dia 
!  Ct=0.01 


!  Weight  1  m  x  220  mm  dia 
!  Ct=0.01 


!  Surface  Anchor 
'  Top  of  Can 

!  Bottom  of  Can 
!  D2 

!  HI 
!  H25 

!  H2  to  H24 
!  D2 

!  Top  of  Weight 
!  Bottom  of  Weight 


1  Tube 


!  Strength  Member  #1 
!  Strength  Member  #2 


!  Tube 


!  Strength  Member  #1 
!  Strength  Member  #2 
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!  Hydrophone  Section  -  HI  to  H25 

TYPE, 3  !  Tube 

MAT,  3 
REAL,  3, 

E,  1 , 2 
TYPE,  2 
MAT,  2 
REAL,  2 

E,i,2  !  Strength  Member  #1 

E,i,2  !  Strength  Member  #2 

EGEN, 24 , 1, -3 


I  * 

!  Lower  VLA  -  H25  to  D2 
TYPE, 3 
MAT,  3 
REAL, 3, 

E, 25,201 
TYPE, 2 
MAT,  2 
REAL,  2 
E, 25, 201 
E, 25, 201 

!  * 

!  Lower  VLA  -  D2  to  Weight 
TYPE, 3 
MAT,  3 
REAL,  3, 

E, 201, 202 
TYPE, 2 
MAT,  2 
REAL,  2 

E, 201, 202 
E, 201, 202 

!  * 

!  Weight 
TYPE, 5 
MAT,  5 
REAL, 5 

E, 202,203 

!  * 

!  * 

Pressure=101 . 3E3 
ESEL, S , MAT,  ,3 
CM , TUBE , ELEM 

SFE, ALL, 1, PRES,  , Pres sure *4 
ESEL, all 

1  Apply  Constraints 

D, ALL,  UY, 0 

D,  101 ,  UX,0,  ,  ,  ,UZ 


!  Tube 


!  Strength  Member  #1 
l  Strength  Member  #2 


!  Tube 


!  Strength  Member  #1 
!  Strength  Member  #2 


!  Atmospheric  Pressure  [PA] 

!  Select  Array  Elements  Material  2 
!  Form  element  component  TUBE 
!  Apply  Internal  Pressure  To  TUBE  Elements 
!  Reselect  all  elements 

l  Fix  all  nodes  in  UY 
!  Fix  anchor  node  101  in  UX  and  UZ 


i  * 


!  Load  with  Gravity 
ACEL, 0,0, 9.81 


i  * 


!  * 

/DSCALE, 1, 1  !  Deflection  Scaling  -  True  Geometry 

/VIEW, 1 ,  , -1 

/PNUM, NODE , 1 

/NUMBER, 0 

/PBC ,  U,  ,1 

/PBC,  ROT,  ,1 

/ PBC , ACEL ,  ,1 

EPLO 

SAVE 

i  * 
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!  Transient  Effects  OFF 
!  Ramped  Loads 


I  Apply  Tension  to  weight 
SS  0,  %tension%  N,  %current%  knots 
!  Solve  Loadstep  #1 


!  0.6  knots 
!  For  Mooring  Rope 


»  For  Tube 
!  For  Can 
!  For  Weight 


I  Turn  ON  Structural  Inertia  Effects 
!  Minimum  Timestep  0.005  s 
!  Maximum  of  100  Equilibrium  Operations 
!  Results  at  end  of  Loadstep 
!  TOLER  =  1%,  MINREF  =  1  N 


!  Increment  time 
!  Apply  tension 


B.2  Solution  Input  File 

!  spread2 .sol 01 
/SOLTJ 

!  Model  Designation  spreadl 
!  parameters 
tension=1000 
dtime=100 
current =0 

i  * 

ANTYPE ,  4 
TRNOPT ,  FULL 
LUMPM,  0 

EQSLV, FRONT ,16-08,0 

AUTOTS , OFF 

NLGEOM , ON 

SSTIF, ON 

NROPT , AUTO 

TIMINT, OFF 

KBC,  0 

j  ★ 

OUTPR,  BASIC,  LAST 
OUTRES , ALL , ALL 
1  * 

TIME,  dtime 
NSUBST, 1 
NCNV,  2 

F, 203, FZ, -tension 

/TITLE, Suspended  VLA  Analysis  - 

SOLVE 

i  * 

I  * 

current  =  0.6 
!  * 

TB, WATER, 1 

TBMODIF, 2,1, -depth 
TBMODIF, 2,2, current *0 . 5144 
TBCOPY , WATER, 1 , 3 
TBCOPY , WATER ,1,4 
TBCOPY, WATER, 1 , 5 

t  * 

AUTOTS , ON 
LNSRCH , ON 
PRED,ON,  , ON 
TIMINT,  ON,  STRUC 
DELTIM, 0.005, 0,0,0 
NEQIT, 100 
OUTRES , ALL, LAST 
CNVTOL , F ,  ,0.01,2,1 

i  * 

*DO, I, 1,10,1 

dt ime =dt ime +100 
F, 203 , FZ, -tension 
TIME, dtime 

/TITLE, Suspended  VLA  Analysis 
SOLVE 

tension=tension-100 

*ENDDO 

FINISH 


-  SS  0,  % tens ion %  N,  %current%  knots 
!  Solve  Loadstep 
I  Decrement  tension 
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